Difference of optical conductivity between one- and two-dimensional doped nickelates 
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We study the optical conductivity in doped nickelates, and find the dramatic difference of the spec- 
trum in the gap (o><4 eV) between one- (ID) and two-dimensional (2D) nickelates. The difference 
is shown to be caused by the dependence of hopping integral on dimensionality. The theoretical 
results explain consistently the experimental data in ID and 2D nickelates, Y2- a; Ca :E BaNi05 and 
La2- a; Sr I Ni04, respectively. The relation between the spectrum in the X-ray aborption experiments 
and the optical conductivity in L^-xSr^NiCU is discussed. 



PACS numbers: 78.20.Bh, 71.10-Fd, 71.20.Be 

Nickelates have received special attention since the dis- 
covery of high T c cuprates. An antiferromagnetic insu- 
lator La2NiC>4 has the same layered structure as that of 
one of the parent compounds of high T c superconduc- 
tors, La2Cu04. In contrast with the cuprate with spin- 
1/2, the nickelate is a spin-1 system. Therefore, the hole- 
doped nickelate, L^-zSr^NiO^ provides an opportunity 
to examine a role of the spin background in the carrier 
dynamics in two-dimension. Another doped nickelate,, 
Y2- 2; Ca 2 ;BaNi05, has also been studied extensivelyBtl 
In the undoped one, Y2BaNiOs, apex-linked NiC>6 octa- 
hedra are arranged in one-dimension and exhibit a spin-1 
chain. 

The spectrum of the optical conductivity in the low 
energy region (w^4 eV) reflects the key parameters 
for the electronic structure such as Coulomb interac- 
tion, Hund's-rule coupling, superexchange interaction, 
hopping integral of carriers, electron-phonon interaction 
and so on. The experiments in one-dimensional (ID) 
nickelate, Y2_ ;E Ca 2: BaNi05, and two-dimensional (2D) 
one, Lap—zSr^NiC^, have been performed by several 
groups .ErQ Both the undoped nickelates are insulators 
with the charge excitation gap of about 4 eV. Upon dop- 
ing of holes, the spectrum appears in the gap. There are 
dramatic differences in the spectrum between ID and 2D 
nickelates. In 2D, the spectrum spreads in the gap (w<4 
eV) and is broader than that in ID. At around 1 eV, the 
spectrum shows two-peak structure in 2D, whereas it has 
a single peak in ID. 

In this paper, we propose that such a remarkable differ- 
ence between ID and 2D nickelates is caused by the hop- 
ping integral of 3d electrons. The hopping integral exists 
between all of the neighboring e g orbitals in 2D. In ID, 
on the other hand, the integral is restricted within the 
3d(3z 2 — r 2 ) orbitals, where the chain direction is taken 
to be the z-axis. It is shown that the dependence of the 
hopping integral on dimensionality affects the dynamics 
of carriers and the spectrum in optical conductivity. We 
examine the optical conductivity in ID and 2D nickelates 
by using the numerically exact diagonalization method. 
We find the dramatic difference of the spectrum between 
ID and 2D: The spectrum in 2D spreads in the gap in 
contrast with that in ID. In 2D, the two-peak structure 
is obtained in the low energy region. One of the peaks is 



due to a low-spin state induced by a doped carrier and is 
governed by the Hund's-rule coupling. The other origi- 
nates in the antiferromagnetic spin background disturbed 
by the carriers. In ID, the peak for the low spin state 
is strongly suppressed. The difference of the spectrum 
is well described by the dependence of hopping integral 
on dimensionality. We also examine the relation between 
the optical conductivity and the single-particle excitation 
spectrum. The experimental data of the optical conduc- 
tivity and X-ray absorption spectroscopy are discussed 
in the light of theoretical results. 

A Ni 2+ ion in a NiOg octahedron has two electrons in 
the e g orbitals, 3d(x 2 — y 2 ) and 3rf(3z 2 — r 2 ). Our starting 
model consists of the following terms; hopping integrals 
between the e g orbitals on nearest-neighbor sites, on-site 
Coulomb interactions and Hund's-rule coupling between 
electrons in the e g orbitals. The Hamiltonian is written 
as, 

H = H t + H 1 . (1) 
The kinetic energy of 3d electrons H t is written as, 

Ht= E tZd^n,^ + h.c), (2) 

(l,m) ,a,fi,iy 

where cj a is the creation operator of electron with spin 
(7 in fi orbital at site I, and is the hopping integral 
between orbital n at site I and v at to. The integral 
depends crucially on the dimensionality of the system 
and is expressed as, 
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in the xy plane in 2D and, 
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in ID, where the chain direction is taken to be z-axis. 
Here, t is the integral between neighboring Zd(x 2 ~y 2 ) 
orbitals and the +(— ) sign in the off-diagonal terms in 
Eq. (|3|) corresponds to the integral in the x(y) direction in 
2D. In 2D, electrons in both e g orbitals contribute to the 
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transport. In ID, on the other hand, an electron in the 
3d(x 2 —y 2 ) orbital is localized in each Ni ion since only 
the integral between neighboring 3g?(3z 2 — r 2 ) orbitals is 
finite. As will be shown below, the dependence of the 
hopping integral on dimensionality causes a crucial role 
in the electronic structure. 

The electron-electron interaction Hi is given by, 

Hi = U y^n;, M ,fn;,^,j. + U' ^ n;, M n; |t) 
i.fi i 

- 2-ftT^(Sz,„ • Si tV + -ni. u n^ v ) 
i 

+ K ^( c lu,A,u,l Cl ^l c W + h - c -), (5) 



where u and v denote the 3g?(3z 2 — r 2 ) and 3d(x 2 — y 2 ) or- 
bitals, respectively, »ty jA , j0 .=cJ j/ljtr cj jAt>0 ., nz iM =n; i([liT +n/ )Aa , 
£7 and J7' are the intra- and inter-orbital Coulomb 
interactions, and K is the Hund's-rule coupling with 
U=U'+2K. 

The O 2p orbitals are not taken into account explic- 
itly, since we are interested in the dependence of optical 
conductivity on dimensionality. When the effect of the 
O 2p orbitals is estimated in the NiOe cluster model, we 
find that the atomic parameters are reduced by ~3Q%, 
From the observed atomic parameters in La2Ni04,u£il 
the values of U and K are taken to be ~5 and ^0.7 
eV, respectively, and the value of t is estimated to be 
t~T 2 d / A~10 _1 eV. Since the observed crystal field split- 
ting of e g orbitals mi^NiC^ is much smaller than the 
Coulomb interaction)!^ we neglect it in this study. The 
importance of phonons jn-.the low energy region (w<0.5 
eV) has been discussed.ET0 Since we are interested in the 
higher energy region, the electron-phonon interaction is 
neglected in the present study. 

We calculate optical conductivity and single-particle 
excitation spectra in the numerically exact diagonaliza- 
tion method on 8 site clusters in ID and y/Exy/S site 
clusters in 2D with periodic boundary condition. The 
values of U/t and K/ 1 used in the calculation are 14 and 
2, respectively. The optical conductivity is written as, 
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where |0) and \n) are the ground state and n-th excited 
state with energies Eq and E™ , respectively, D denotes 
the Drude constant, and N$ and N are the system size 
and number of electrons, respectively. The Drude weight 
is estimated from the sum rule l 12 H 13 l The a-axis compo- 
nent of the current operator j a is given by, 
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FIG. 1. Optical conductivity a(ui) in ID and 2D. (a), 
(d): non-hole doping, (b), (e): two hole (25% hole) doping, 
and (c), (f): four hole (50% hole) doping. The parameters 
U/t— 14, K/t=2 are used. The 5-functions are convoluted 
with a Lorentzian broadening of 0.2t. The Drude constants 
(D) are not displayed. The values of D are; (a):0.00, (b):2.96, 
(c):6.68, (d):-1.27, (e):3.10, and (f):2.50 in the same unit. 



where a—z and x in ID and 2D, respectively. In the same 
way, the single particle excitation spectrum is expressed 
as, 

A,(k,w) = |(n|c fe) ^|0)| 2 5(-a; - E^ 1 + <) 

71 

+ ElH4,^|0)| 2 5(a;-< +1 +O, (8) 

n 

where Cfc iM a is the Fourier transform of ci 4l-a . 

In Fig. |l|, the numerical results of cr(u>) in ID and 2D 
with no- hole (insulating), two holes (25% doping) and 
four holes (50% doping) are shown. In the undoped 
states, a charge excitation gap exists at w/i~12. When 
the value of t is taken to be ~0.33 eV, the gap energy 
is estimated to be^A eV, which is consistent with the 
experimental dataira Uppa doping of holes, the spec- 
trum appears in the gapJij The shape of the spectrum 
is quite different between ID and 2D systems. In 2D, 
the spectrum has some characteristic peaks in the gap 
(w/t<12), whereas in ID it is rather sharp in the low 
energy region(w/t~2). 

To clarify the nature of the spectrum in the gap, let 
us first consider the charge excitations in the strong cou- 
pling limit (i<C U, K). In the undoped states, each site 
has two electrons in the e g orbitals with high spin (5=1) 
configuration due to the Hund's-rule coupling. There- 
fore, the charge excitation gap has the energy ~L7. In 
the hole-doped systems, there exist singly occupied sites 
with 5=1/2. The hopping of an electron from a dou- 
bly occupied (5=1) site to a singly occupied (5=1/2) 
site costs the energy less than the gap energy. After the 
hopping, the high spin (5=1) or low spin (5=0) state 
appears at the doubly occupied site. In the former case, 
the excitation is caused by the hopping of electron which 
disturbs the spin background. In the latter case, on the 
other hand, the excitation is characterized by the energy 
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FIG. 2. Single particle excitation spectra A l _ t (k,ui) with 
four holes (50% hole doping case). Left and right panels are 
for ID and 2D, respectively. Upper and lower panels are 
the spectra for electrons in 3d(3z — r ) and 3d(x —y ) or- 
bitals, respectively. The dotted lines denote Fermi level. The 
<5-functions are convoluted with a Lorentzian broadening of 
0.4t. The parameter values are the same with those in Fig. 1. 



scale with K. These excitations appear in the gap of 
a(uj). 

Let us next calculate the single particle excitation 
spectrum A^(k,cu) and examine the relation between 
A^k, oS) and a(u>). Figure |2| shows A^{k,uj) for 50% 
hole doped systems in ID and 2D. The peaks (i)~(iii) 
are characterized by the 3 ' (high spin), l E (low spin), 
and 1 A\ (low spin) states at the site where an electron 
is added, respectively, and the peak (iv) is given by the 
triply occupied states at the site. In Fig. g(a), the peak 
(iii) does not exist, because the 1 A± state, in which two 
electrons are in the same e g orbital, can noJ_.be created by 
adding an electron in 3d(3z 2 —r 2 ) orbital.E__l In Fig. ||(b), 
the spectrum is independent of k, since the hopping ma- 
trix elements for 3d(x 2 — y 2 ) orbitals do not exist in ID 
and electrons in the orbitals are localized in Ni ions. In 
Figs. ||(c) and (d), the hopping matrix elements for both 
orbitals cause the broadening of the spectra of (ii)~(iv). 

In ID, since only 3d(3z 2 — r 2 ) orbital contributes to the 
current, the excitations from the states below the Fermi 
level to those of (i) and (ii) in Fig. ||(a) appear in the 
gap of a(uo). On the other hand, in 2D, the excitations 
to the states of (i), (ii), and (iii) in both Figs, ^(c) and 
(d) appear in the gap, so that the spectrum is broader 
in 2D than that in ID. To resolve these excitaions in the 
optical conductivity, we introduce the operator j a given 
by, 

ja = —i t'i +a jPLacj flcr (ci-a^,a—Cl+a,iy,a), (9) 

where Pi tCC is the projection operator which restricts the 
state to a at site I. For example, the projection operator 
for a^Ai with two electrons is written as, 
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FIG. 3. Optical conductivity are analyzed by using the 
operators given by Eq. Left and right panels are for ID 
and 2D, respectively, (a), (c): two hole (25% hole) doping, 
(b), (d): four hole (50% hole) doping. The thin lines are 
the same results shown in Fig. The dot-dashed, thick, 
gray, and dashed lines denote spectra for a— 3 A2, 1 E, l Ai, 
and the triply occupied state, respectively. The 5-functions 
are convoluted with a Lorentzian broadening of 0.2t. The 
parameter values are the same with those in Fig. 1. 



+ ?t*,t>,T n 2,ta(i - ™z,«,t)(i - n i,u,i) 

+ c J t ,u,T C /,u,i Cj .«.l Ci .".t + C t,v,1 c l,v,l c l,u,lCl,u,l}- ( 10 ) 

We calculate the spectrum Eq. (0) by replacing j a with 
j a . The results for a— 3 A 2 , 1 E, A\, and triply occupied 
state are shown in Fig. 3. The dot-dashed, thick, gray, 
and dashed lines denote the spectra for a— 3 A2, 1 E, 1 Ai > 
and the triply occupied state, respectively. The thin lines 
show the same results as those in Fig. [l] As shown in 
Fig. 3, the peaks (i), (ii), and (iii) are well described by 
the final states with a= 3 A2, 1 E, and *A\, respectively. 

The spectra of (i) in both ID and 2D are brought about 
by the motion of carriers which disturbs the spin config- 
uration. They are similar to the mid-gap spectra in high 
T c cuprates. Note, however, that the spectra exist even 
in ID 5=1 systems. This is in contrast with thai, in ID 
5=1/2 systems where such spectra do not exist.Ej 

The spectra of (ii) and (iii) are due to the excitations 
in which the motion of carriers changes doubly occupied 
sites from high spin to low spin states. These spectra 
in ID are smaller in magnitude than those in 2D. In 
particular, the spectrum of (iii) does not exist in ID, 
because the hopping matrix element for ?>d(x 2 ~ y 2 ) or- 
bitals does not exist and the 1 A\ state is excluded in 
the optical process. The small magnitude for the spec- 
trum of (ii) in ID is understood as follows; since the hop- 
ping matrix element is restricted within the 3d(3z 2 — r 2 ) 
orbitals in ID, the ferromagnetic spin alignment is sta- 
bilized around a singly occupied site due to the double 
exchange interaction (see Fig. §(a)). In 2D, on the other 
hand, the off-diagonal hopping matrix elements induce 
the antiferromagnetic superexchange interaction and re- 
duce the stability of ferromagnetic spin alignment (see 
Fig. |](b)). Therefore, the probability for the excitation 
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FIG. 4. The difference of the hopping matrix elements 
between ID and 2D. (a) In ID, the hopping matrix ele- 
ments are restricted within 3d(3z 2 — r 2 ) orbitals. (b) In 2D, 
the off-diagonal hopping matrix elements induce the low spin 
state. 



to the low spin state is smaller in ID than that in 2D, 
and the intensity of the peak (ii) is reduced in ID. 

Experimentally, it has been shown that the spetrum 
of the optical conductivity in La2- x Sr x Ni04 spreads in 
the gap and has two peaks at 0.6 and 1.5 eV.H The two 
peaks grow with doping. Figs. 3(c) and (d) explain this 
fact; two peaks (i) and (ii) grow with doping, and the 
energies are ~2i and ~5i which are estimated to tae-r^O.6 
and ~1.5 eV, respectively, where i~0.3 eV is takcnila The 
appearance of the spectrum of (hi) is also consistent with 
the observation that the spectrum spreads in the gap. In 
ID, the peak (i) in Fig. 3(a) is much larger than (ii), 
and the peak (i) explains well the single peak structure 
in Y 2 _ x Ca x BaNi0 5 .El 

Let us discuss the O Is near-edge X-ray absorption 
spectroscopy (XAS). Pellegrin et alS3 have observed that 
there appear the peaks at 528.7 and 530 eV in the po- 
larized XAS upon doping of holes. They identified the 
two peaks as the high spin state and the low spin one, 
respectively. In our theory, the states correspond to the 
states of (i) and (ii) in A v (k, a;). EH It is shown that the 
excitations (i)~(iv) in A fJ _(k 1 Lu) bring about the peaks 
(i)~(iv) in (j(oj), respectively. Thus, we find that the 
peaks at 0.6 eV and 1.5 eV in the optical conductivity 
in La2- a; Sr 2: Ni04 correspond to the excitations at 528.7 
and 530 eV in XAS, respectively. 

In summary, we have examined the optical conductiv- 
ity in ID and 2D nickelates and found that the dramatic 
difference occurs between them, since the hopping in- 
tegral of 3d electrons between neighboring Ni sites has 
the off-diagonal term in the e g orbitals in 2D nickelates, 
whereas the integral is restricted within the 3rf(3z 2 — r 2 ) 
orbitals in ID ones. The theoretical results explain con- 
sistently the experimental data in Y2_ a; Ca a ;BaNi05 and 
La2_ 2; Sr 2 ;Ni04. We have also shown the relation be- 
tween the spectra in XAS and optical conductivity in 
L^-zSr^NiC^. 
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